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SUMMARY 

Contact-impact proBlems involving finite deformation axi symmetric membranes 
are solved By tfie finite element raetbod with explicit time integration. The 
formulation of the membrane element and the contact constraint conditions are 
discussed in this paper. The hyperelastic, compressible Blatz and Ko material 
is used to model the material properties of the membrane. Two example problems 
are presented. 


INTRODUCTION 

The purpose of this paper is to present a method for the dynamic analysis 
of contact-impact problems involving hyperelastic compressible membranes. A 
strain energy functional developed by Blatz and Ko (ref, 1) is used to charac- 
terize the material of the raemBrane. This element was added to HONDO (ref. 2), 
a finite element code that explicitly integrates the equations of motion. The 
contact-impact algorithm, which was also added to HONDO, was recently developed 
by Hallquist (ref. 3) and is briefly described here. 

Two examples are provided to demonstrate the capability of the method; in 
the first, a flat circular membrane is inflated by a pressure loading into a 
thick-walled sphere*, and in the second, the sphere is impacted into the mem- 
brane. 


FORMULATION 
Equation of Motion 

Since an explicit time integration scheme is being considered, the equation 
of motion becomes 

= ^ - f ( 1 ) 

where is the diagonal (lumped) mass matrix^tf is a global vector of nodal ac- 
celerations, P is the applied load vector, and F is the stress divergence vector. 
This equation is integrated by the velocity-centered central difference method. 


Work was performed under the auspices of the United States Energy Research 
and Development Administration under contract No. W-74-05-eng-48. 
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Material Properties 

The strain energy density per unit undeformed volume Ug for a compressible 
hyperelastic material is expressed as 


u, = p 


3 + 



( 2 ) 


where y is the shear modulus, v is Poisson's ratio, and is the ith strain 
invariant. These invariants can be expressed in terms of the principal stretch 
ratios X-| , X 2 , X 3 in the meridional , circumferential , and transverse directions * 
respectively, as 


2 2 2 
+ aI + A3 

2 2 2 
I3 = Af A^ A^ 


(3) 


For thin membranes, the stress component normal to the midsurface is assumed to 
be zero; hence, A 3 can be expressed in terms of A] and A 2 


V 


^3 CA1A2) 


and the strain energy density becomes a function of A] and A 2 . 


(4) 


Membrane Element 

An isoparametric axi symmetric membrane element is shown in Figure 1. The 
R, Z, and meridional coordinates S of the undeformed configuration are related 
to the natural coordi nate L through 


R ^ 1 (1 - L)R^’ + I (1 + L)rJ 

Z = i- 0 - L)Z^ + |- (1 + L)Z^ (5) 

S = ^ (1 - L)sV+ ^ (1 + L)S^ 


and similarly for the displacement components Uy, and Uj 

u^ = i (1 - L)u^ + (1 + L)uj 

In the deformed configuration, the r and z coordinates along the midsurface are 
given by . 


r = u^ + R 
z = u^ + Z 

The principal sketch ratios A-| and A 2 can be defined as 


(7) 
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( 8 ) 



Substitution of equattons (5) and (6) into equation (7), putting the result 
into equation (8), and applying the chain rule leads to expressions for X-| 
and X 2 in terms of the nodal point quantities 



(1 ^ Dul + 0 - DuJ 

X = 1 + j ■ — — — ■ ■ • 

(1 - L)R^ + (1 + L)r 


where t = SJ S'*'. 


Since X-j and X 2 are now functions of the natural coordinate L, the total 
strain energy stored within the membrane element during deformation can be ex- 
pressed as the integral 



( 10 ) 


in which h is the thickness of the undeformed membrane. 

The partial derivatives of U with respect to the nodal displacement compo- 
nents yield nodal point forces that are subsequently accumulated into the stress 
divergence vector. In the problem aqder consideration these derivatives can be 
calculated very easily. For example, the nodal point force acting in the r- 
direction at the ith node is given by 



where T] and T 2 are Lagrange stresses in the meridional and circumferential di- 
rections, respectively, A two point Gauss quadrature is used to perform the 
above integrations. 


The lumped masses for each element are found by the addition of the off- 
diagonal terms of the consistent mass matrix to the diagonal term. Each mem- 
brane element yields the following contributions to the nodal point mass at 
nodes i and j, respectively. 


m^. = 2wpXh (rVs + R^/6) 


( 12 ) 


m. = 2TTpilh (rJ’/ 3 + RVe) 

vl 


where pis the mass density of the undeformed membrane. 


For stability the time step At is restricted such that the inequality 
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2 4 

At < -j (13) 

is satisfied where X^ is the maximum eigenvalue of in which K is the stiff- 
ness matrix. 


A time step At is calculated for every element in the mesh and 90 percent 
of the smallest value is then used. For the membrane element X^ is calculated 
exactly from 





2 

9^U 


"a 


,3u 


J2 



(14) 


Contact-Constraint Condi tions 

Two elastic bodies occupying regions b 1 and b 2 in the reference configura- 
tion at time t = 0 are shown In h'gure 2. The boundaries of b 1 and B^ are de- 
noted by 3$! and respectlyely. After deformation at time t ^ 0, these 
bodies occupy regions bl^and b^A The boundaries of b^ and b^ are denoted by 
3sl and 3s2. Whenever bl and b^ are in contact, the nodal points on 3$! in the 
contact region are constrained to slide along line segments connected by nodal 
points lying on 3s2. Separation Is permitted when the interface pressure Is 
negative. Impact and release conditions are applied whenever nodal points on 
3s* come into contact with 3s2. These conditions, which are based on the gen- 
eralization of those given by Hughes, et al . (ref. 4). conserve linear and angu- 
lar momentum. 


Constraint conditions must be Imposed into the equations of motion for each 
node of 3sl in contact with a segment of 3s^. These conditions are imposed 
through a transformation of displacements which is performed at the beginning of 
each time step. In this transformation the radial and vertical displacement com- 
ponents of the node on 3sl are transformed into a displacement component tangen- 
tial to the segment and a relative displacement component normal to the segment 
of 3s^ on which It rests. Since no separation is permitted during the time step 
the displacement, velocity, and acceleration of this latter component are set to 
zero. A transformation matrix T is constructed which relates the vector of 
global displacements l!f to a vector i!f' containing the transformed components 

ti = f ti' (15) 


Letting t remain constant throughout the time step and differentiating equation 
(15) with respect to time yields 


T- 

u - T u 


(16) 


Equation (16) is substituted into equation (1) and the resulting equation is 
premultiplied by T^ In order to obtain the modified equations of motion 

i^' = f*(f - f) (17) 

whlch^contains the contact constraints. Here M* = T^ iff. Although is diag- 
onal m' is not. For computational efficiency the appropriate off-diagonal mass- 
es are summed to the diagonal. 

After eguatlon (17) is solved for iJ'', the normal accelerations of the nodes 
of 3 si on 3 s 2 relative to 3s^ are set to zero. The global accelerations then 
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follow directly from equation (16), 


EXAMPLES 

In the following examples, all physical quanttttes are given in nondimen- 
sional form. Any consistent units may Be assumed without altering the results. 

Inflation of a Membrane into a Thick’-Walled Sphere 

A flat unstretched circular membrane with a thickness of 0,01 and a radius 
of 2.0 is positioned beneath a thick-walled sphere having an inner radius of 
0.40 and an outer radius of 0,60, In the undeformed configuration, the distance 
measured perpendicularly from the center of the monbrane to the center, of the 
sphere is 1.20, The hyperelastic material described by equation (2) is used to 
model the material of both the membrane and the sphere with v and -p set to 
0.463 and 150, Densities of 1.5 and 0.15 were assumed for the material of the 
membrane and sphere, respectively. 

The membrane is inflated By a pressure p defined by 
0 < t < 0,11 p = 1.250 

0.11 < t < 0.15 p = 1.250 - 1.125 (18) 

t > 0.15 p = 0.125 

and is brought into contact with the sphere. 

In Figure 3 the deformed shapes at various times throughout the deforma- 
tion time history are shown. At late times some wrinkling occurs (for example, 
note the last frame) and the calculation ceases to be meaningful, A total of 
eighty elements were used in the calculation. Forty elements were of the mem- 
brane type. 


Thick-Walled Sphere Impacting a Membrane 

In this example the thick-walled sphere impacts the flat circular membrane 
with an initial velocity of 1.0. The dimensions and material properties of the 
membrane and sphere are identical to those of the preceding example. In Figure 
4 the deformed shapes at various times are shown. Maximum deflection occurs at 
the center of the membrane at approximately t = 0,90 after which rebound begins. 
Separation of the sphere and membrane occurs at approximately t = 1.94. 

In the above examples the stress at the center of the membrane increases 
significantly after the initial contact thereby providing -evidence that a large 
amount of sliding occurs during contact. 
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t = 0.0 


t = 0.045 


1=0088 
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t = 0.233 t = 0.264 1=0.296 

Figure 3.- Inflation of circular membrane into 
thick-walled sphere* 
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Figure 4.- Impact of thick-walled sphere into 
circular membrane. 



